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Time  series  of  bio-optical  properties  in  a  subtropical  gyre: 

Implications  for  the  evaluation  of  interannual  trends 
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[i]  With  a  validated  Quasi-Analytical  Algorithm,  an  1 1  year  (1998-2008)  monthly  time 
series  of  the  primary  optical  properties  of  waters  in  the  center  of  the  South  Pacific  gyre  was 
developed  from  Sea-viewing  Wide  Field-of-view  Sensor  (SeaWiFS)  and  Moderate 
Resolution  Imaging  Spectroradiometer  (MODIS).  Also  derived  are  chlorophyll  a  (Chi  a) 
concentrations  with  the  operational  empirical  algorithms  for  SeaWiFS  and  MODIS  The 
optical  properties  include  the  absorption  coefficient  (at  443  nm)  of  phytoplankton  (aph)  and 
that  of  the  combination  of  detritus  and  gelbstofT  (adg).  From  these  time  series,  we  further 
derived  their  annual  background  (summer  low)  and  seasonal  intensity  (the  difference 
between  winter  high  and  summer  low).  These  time  series  show  that  ( 1 )  the  optical  properties 
have  different  seasonal  and  interannual  variations,  indicating  different  dynamics  of  these 
properties  in  the  subtropical  gyre;  (2)  there  is  a  decreasing  trend  (r2  =  0.24)  of  the  background 
ap h  in  the  1 998  2008  period  and  an  increasing  trend  of  the  aph  seasonal  intensity  (r  0  11) 
for  this  period,  and  both  trends  are  not  statistically  significant;  (3)  the  aph  time  series 
agrees  with  the  Chi  a  time  series  at  the  seasonal  scale,  but  differs  with  respect  to  interannual 
variations;  and  (4)  different  interannual  trends  could  be  inferred  with  different  time 
frames.  These  results  emphasize  that  it  is  difficult  to  draw  unequivocal  conclusions  about 
long-term  trends  of  biogeochemical  properties  in  the  oceans  with  the  current  relatively  short 
bio-optical  records.  To  clarify  and  predict  such  trends,  it  is  critical  to  get  a  full  account  of  the 
forces  that  are  responsible  for  the  seasonal  and  interannual  variations  of  these  properties. 
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in  a  subtropical  gyre:  Implications  for  the  evaluation  of  interannual  trends  of  biogcochcmical  properties,  J  Geophys.  Res. ,  115, 
C090 1 2,  doi:  1 0. 1 029/2009JC005865. 


1.  Introduction 

[2]  Subtropical  gyres  arc  not  only  important  for  studies  of 
ocean  biogcochcmistry  and  the  global  carbon  cycle 
[Behrenfeld  et  al. ,  2006;  McClain  et  al.,  2004],  they  arc  also 
important  for  the  detection  of  anthropogenic  climate  change 
[Barnett  et  al .,  2001  ].  Recent  studies  based  on  chlorophyll  a 
(Chi  a)  concentration,  derived  from  satellite  ocean  color 
measurements,  have  shown  that  Chi  a  in  the  subtropical  gyres 
is  showing  a  waning  trend  [Gregg  et  ai ,  2005]  and  that  the 
areas  of  low'  Chi  a  gyres  are  expanding  [McClain  eta  l ,,  2004; 
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Polovina  et  a l .,  2008].  These  results  [Gregg  et  al .,  2005; 
Polovina  et  al. ,  2008]  imply  a  reduced  production  of  phyto¬ 
plankton  in  surface  waters  of  the  gyres  in  the  past  decade  and 
contribute  to  speculation  that  such  an  annual  declining  trend 
could  be  a  result  of  enhanced  surface  water  stratification  due 
to  temperature  increases  [Gregg  et  al .,  2005,  Polovina  et  al ., 
2008]  However,  each  of  these  studies  used  a  single  time 
window  to  estimate  long  term  trends  with  the  consequent 
possibility  of  alias.  Further  when  analyzing  satellite  ob¬ 
servations,  these  studies  used  Chi  a  that  was  derived  with  an 
empirical  algorithm  developed  and  validated  with  global  data 
sets  [O  'Reilly  et  al. ,  1998]  and  not  specifically  validated  for 
waters  in  the  gyres.  Such  empmeal  algorithms  use  the  blue 
green  ratio  of  ocean  color  to  derive  Chi  a ,  where  the  poten¬ 
tially  confounding  impacts  of  gelbstofT  (or  colored  dissolved 
organic  matter  (CDOM))  could  not  be  removed  [Carder 
et  al.,  1991]. 

[3]  Here,  we  use  in  situ  bio-optical  data  measured  near  the 
surface  center  of  the  South  Pacific  gyre  to  show  that 
absorption  coefficients  of  gelbstofT  and  phytoplankton  of 
sueh  extremely  clear  waters  can  be  derived  with  high  accu¬ 
racy  from  remote  sensing  reflectance  with  a  scmianalytical 
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Figure  1.  Location  of  the  study  site  (23-28°S,  130-  140°W). 
The  background  image  is  the  annual  water  clanty  (measured 
by  cuphotic  depth  in  meters;  see  the  work  of  Lee  et  ai  [2007] 
algorithm)  for  2008  (from  Ocean  Biology  Processing  Group). 

algorithm.  We  then  use  this  algorithm  to  develop  an  1 1  year 
(January  1 998-Deccmber  2008)  monthly  time  series  of 
absorption  cocfTicicnts  for  phytoplankton  and  gclbstoff  based 
on  ocean  color  radiometric  data  from  the  Sea  viewing  Wide 
Ficld-of-view  Sensor  (SeaWiFS)  and  the  Moderate  Resolu¬ 
tion  Imaging  Spectrometer  (MOD1S)  (both  operated  by 
NASA  [McClain,  2000]).  A  monthly  time  scries  of  Chi  a  was 
also  derived  from  the  same  ScaWiFS  and  MODIS  data  with 
the  operational  empirical  algorithms  (OC4v4  and  OC3m) 
[O  'Reilly  et  al ,  2000].  Seasonal  and  interannual  variations 
were  evaluated,  and  mterannual  trends  were  analyzed  with 
different  time  frames.  Results  from  these  studies  indicate  that 
(1)  a  long-term  (decades)  record  of  ocean  color  observations 
is  required  to  deduce  unequivocal  trends  of  biogcoehcmical 
properties  in  the  oceans  and  (2)  it  is  more  important  to 
account  for  the  forces  that  arc  responsible  for  the  seasonal 
and  intcrannual  variations,  as  such  forces  arc  the  foundation 
to  understand  and  eventually  predict  the  apparent  trends, 
if  any. 

2.  Data  and  Method 

2.1.  In  Situ  Data  for  Algorithm  Validation 

[4]  Remote  sensing  spectral  reflectance  (/?re(A);  sr  *;  ratio 
of  water  leaving  radiance  to  downwclling  irradiancc  above 


sea  surface)  for  five  stations  (STB4,  STB5,  STB8,  GYR4,  and 
GYR5  that  arc  listed  in  Table  1  of  Morel  et  al  [2007])  in 
the  South  Pacific  gyre  was  computed  from  direct  mea¬ 
surements  of  downwclling  irradiancc  above  the  sea  sur¬ 
face  (£5(A);  W  m  nm  *)  and  nadir  upwclling  radiance 
(LU(A);  W  m  nm  'sr  ).  Measurement  of  £*( A)  was  taken 
aboard  ship  while  measurement  of/.M(  A)  was  made  at  a  depth 
of  20  cm  below  the  ocean  surface  using  a  modified  hyper- 
spectral  profiling  radiometer  (HyperoPro,  Satlantic,  Inc.) 
adapted  to  float  at  the  sea  surface  and  tethered  such  that  the 
instrument  operated  at  a  distance  ofHOO  m  from  the  vessel. 
Instrument  tilt  was  measured  directly,  and  measurements 
were  rejected  if  tilts  exceeded  5°.  Measurements  were  made 
over  the  wave  band  380-  800  nm  with  a  spectral  resolution  of 
3.3  nm,  with  each  band  having  a  half  maximum  band  pass 
of  10  nm.  Dark  values  were  taken  every  five  samples  by  use 
of  an  internal  shutter.  These  were  linearly  interpolated  for 
each  light  value,  and  then  subtracted  from  the  observations. 
Calibration  coefficients  and  corrections  for  immersion  effects 
were  obtained  following  standard  protocols  [Mueller  et  al. , 

2003]  and  applied  to  the  measurements;  demonstrated  abso¬ 
lute  accuracies  arc  <2.8%  for  radiance  and  <2.1%  for  irra¬ 
diancc  [Gordon  et  al. ,  2009]  Irradiancc  and  radiance  data 
were  taken  for  3  min  at  each  deployment,  with  each  obser¬ 
vation  within  the  deployment  time  series  representing  inte¬ 
gration  times  of0.03  to  0.5  s,  depending  on  the  intensity  of  the 
incident  (ir)radiancc.  These  measurements  were  then  inter¬ 
polated  to  a  common  time  frame  at  a  frequency  of  every  2  s 
and  to  a  common  spectral  resolution  of  every  2  nm. 

[5]  Upwclling  radiance  measurements  were  then  propa¬ 
gated  to  the  sea  surface  using  an  iterative  approach  [Mueller 
et  al ,  2003]  that  estimates  the  spectral  diffuse  attenuation 
coefficient  from  spectral  ratios  of  measured  radiance,  and  the 
nadir  water-leaving  radiance  above  the  sea  surface,  Z,W(A),  is 
then  computed  based  on  Fresnel  reflectance  and  the  real 
relative  index  of  refraction  [Mueller  et  al. ,  2003].  A  3  min 
time  series  of  R was  made  by  dividing  the  computed  water 
leaving  radiance  by  the  downward  irradiancc  for  each  time 
interval,  and  an  average  value  and  standard  deviation  was 
computed  for  each  deployment. 

[6]  Concurrent  data  of  inherent  optical  properties  (IOP:  the 
total  absorption  coefficient  (a;  m  *)  and  the  portion  attribut¬ 
able  to  particulates  (ap;  m  1 )  at  420  nm)  of  these  stations  were 
measured  and  presented  in  Table  1  in  the  work  of  Morel  et  al 
[2007].  The  absorption  coefficient  attributable  to  yellow 
substance  or  gclbstoff  (ag;  m  ')  at  420  nm  is  estimated  from 
the  presented  ag(310)  value,  with  a  spectral  slope  equal  to 
0.0203  nm  1  (for  the  wavelength  range  of  300  500  nm) 
reported  in  the  same  article.  No  particle  backscattcnng 
coefficients  (h*^)  arc  listed,  but  they  can  be  inferred  from  the 
work  of  Huot  et  ai  [2008]  and  Twardowski  et  al.  [2007]  that 
hb p(550)  was  in  a  range  of  0.0002  0.0004  m  1  based  on  the 
listed  Chi  a  values. 

2.2.  Satellite  Data  for  Temporal  Studies 

[7]  The  research  site  selected  to  study  the  temporal  varia¬ 
tions  ofoptical  and  biological  (Chi  a)  properties  is  a  small  box 
(23-28°S,  130-140°W)  in  the  South  Pacific  gyre  (see 
Figure  1 ),  which  approximates  the  surface  center  of  this  gyre 
where  the  spatial  distribution  of  bio  optical  properties  arc 
nearly  uniform  [Behrenfeld  et  al. ,  2005;  McClain  et  al. , 

2004] .  This  site  was  also  selected  because  (1)  the  waters  arc 
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Figure  2.  Comparison  of  remote  sensing  reflectance  between  SeaWiFS  (solid  squares)  and  MODIS 
(open  circles)  measurements. 


the  ‘‘clearest”  of  the  world’s  oceans  [Morel  et  aL ,  2007],  thus 
it  is  easier  to  detcet  changes  of  biogeoehemical  properties, 
if  any,  and  (2)  the  spatial  gradient  is  minimum  [ Behrenfeld 
el  al. ,  2005],  thus  changes  observed  here  approximate  local 
changes  instead  of  results  of  advection.  In  addition,  there 
are  recent  bio-optieal  measurements  measured  in  this  gyre 
(sec  seetion  2.1)  to  validate  algorithms  for  remote  sensing 
retrievals. 

[x]  For  all  pixels  in  this  box,  monthly  binned  normalized 
water-leaving  spectral  radiance  (Lwn(A))  data  with  a  0  km 
resolution  were  downloaded  from  the  OccanColor  Web 
portal  (http://oeeaneolorgsfc.nasa.gov)  that  is  operated  by 
the  Ocean  Biology  Processing  Group  supported  by  NASA. 
These  Lwn(^)  values  include  those  from  SeaWiFS  (412, 443, 
490,  5 1 0,  555,  and  667  nm;  January  1 998  to  December  2007; 
after  reprocessing  5.2)  and  those  from  MODIS-Aqua  (412, 
443, 488,  53 1 , 55 1 ,  and  667  nm;  from  July  2002  to  December 
2008;  after  reprocessing  1.1).  In  order  to  minimize  impacts  of 
random  spikes  on  the  inverted  products  and  analysis  that 
followed,  for  each  month  (and  each  sensor)  the  top  and  bot¬ 
tom  10%  of  the  values  within  the  box  were  excluded  before  a 
spatially  averaged  Lwn(A)  was  calculated.  We  also  compared 
the  median  value  with  the  average  for  each  band  and  found 
that  the  two  were  nearly  identical  (the  differences  were  under 
0.5%).  Therefore,  this  average  value  was  used  to  represent 
the  property  of  the  study  site  for  that  sensor.  The  multiband 
LWn(A)  values  were  then  converted  to  multiband  remote 


sensing  reflectance  (/?re(A))  by  dividing  LWN(A)  by  the  extra¬ 
terrestrial  solar  spectral  irradianee  (/*o(  A))  [ Gordon ,  2005]. 

[9]  The  remote  sensing  reflectance  values  at  667  nm  from 
SeaWiFS  (/?rs(667))  were  around  0.00026  sr  !,  which  arc 
considerably  higher  than  that  from  MODIS  (~0.00012  sr  *) 
for  such  clear  waters  [Franz ,  2009].  As  Rn(667)  is  generally 
between  0.00005  and  0.00016  sr  1  for  particulate  backscatter 
values,  />bp(667),  between  0  and  0  001  m*"1  based  on 
Hydrolight  [Mobley,  1995]  simulations,  these  high  values 
suggest  a  bias  (overestimation)  of  SeaWiFS  /?re(A)  [also  see 
Franz ,  2009]  for  these  waters  (note  that  bhv( 532)  here  is 
MJ.0002  m  1  [Twardowski  el  al.  2007]).  To  be  more  con¬ 
sistent  with  theoretical  predictions  and  those  from  MODIS, 
we  adjusted  the  SeaWiFS  Rn( A)  by  subtracting  an  offset  of 
0  000175  sr  1  for  all  bands  (0.000175  sr”1  was  deduced  by 
forcing  the  average  of  SeaWiFS  Rrs(66l)  to  be  equal  to  the 
average  of  MODIS  /?rs(667)). 

[10]  Figure  2  presents  the  time  senes  of  adjusted  SeaWiFS 
Rn  at  4 1 2,  443,  55x  (555  for  SeaWiFS,  55 1  for  MODIS),  and 
667  nm  along  with  those  from  MODIS  for  the  2002-2007 
period.  Clearly,  the  two  sets  of  R ^  match  each  other  very  well 
f-2.2%  average  percentage  difference)  at  the  4 1 2  and  443  nm 
bands  (the  two  spectral  channels  that  arc  critical  for  analyti¬ 
cally  separating  contnbutions  of  gclbstoff  and  phytoplank¬ 
ton),  although  slightly  larger  (-5%)  differences  do  exist  at 
412  nm  in  the  Austral  summer  of  2007.  The  consistency  of 

provides  us  a  confidence  to  achieve  longer  time  series  by 
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combining  and  merging  ScaWiFS  and  MODIS  observations. 
ScaWiFS  /?re(5 55)  is  slightly  (-7.6%)  lower  than  that  of 
MODIS  /?re(551)  because  of  the  slightly  different  band  cen¬ 
ters,  and  its  seasonal  pattern  is  subdued.  Compared  to  MODIS 
/?re(667),  ScaWiFS  /?re(667)  is  quite  noisy  with  no  apparent 
seasonal  patterns.  These  reduced  seasonal  features  in  are 
due  to  the  smaller  digitization  ( 1 0  bit)  of  the  ScaWiFS  sensor 
[Hu  et  al.,  2001]. 

[u]  For  analysis  and  comparison,  we  downloaded  the  sea 
surface  temperature  (SST)  products  released  by  the  Physical 
Oceanography  Distributed  Active  Archive  Center  at  the 
NASA  Jet  Propulsion  Laboratory  (Pasadena,  CA,  http:// 
podaac.jpfnasa.gov),  which  was  measured  by  NOAA 
AVHRR  Pathfinder  and  has  a  spatial  resolution  of  0  x  9  km 
A  spatially  averaged  SST,  similar  to  the  processing  of  LWN 
data,  was  obtained  to  represent  the  monthly  time  scries  from 
1998  to  2008. 

2.3.  Inversion  Method 

23.1.  Correction  of  Raman  Scattering 

[12]  Water-leaving  radiance  is  composed  of  contributions 
from  clastic  and  inelastic  scatterings  [Lee  et  al.y  1994; 
Mobley ,  1994;  Sathyendranath  and  Platt,  1997].  The  latter 
part  includes  Raman  scattering  [Stavn  and  Weidemann , 
1988],  and  fluorescence  from  chlorophyll  [Gordon,  1979; 
Yentsch  and  Phinney ,  1985]  and  gelbstoff  [Lee  et  al. ,  1994; 
S.  K.  Hawes  et  al.,  Quantum  fluorescence  efficiencies  of 
marine  humic  and  fulvic  acids:  Effects  on  ocean  color  and 
fluoromctric  detection,  paper  presented  at  Ocean  Optics, 
1992].  In  this  study,  we  ignored  fluorescence  because  both 
were  very  weak  (<I%  to  R „)  for  such  extremely  clear  waters. 
On  the  other  hand,  since  water  constituents  are  scarce  in  the 
gyres  [Morel  et  al 2007],  there  is  sufficient  blue-green 
sunlight  in  the  water  column  to  stimulate  Raman  scattering, 
whose  contribution  to  R is  then  not  negligible  [Hu  and  Voss, 
1997;  Marshall  and  Smith,  1990;  Sathyendranath  and  Platt, 
1998;  Stavn  and  Weidemann,  1988].  This  effect  needs  to 
be  accounted  for  when  an  analytical  inversion  algorithm  is 
applied.  Unlike  the  /?re  models  for  clastic  scattering  (whieh 
arc  generally  functions  of  lOPs)  [Gordon  et  a/,,  1988; 
Sathyendranath  and  Platt,  1997],  Raman  contribution  to  R ^ 
is  not  only  a  function  of  lOPs,  but  also  a  function  of  incident 
irradiancc  at  the  input  and  output  wavelengths  [Lee  et  al , 
1994;  Marshall  and  Smith ,  1990]  For  R from  ScaWiFS 
and  MODIS  measurements,  however,  no  such  data  are 
available.  To  overcome  the  shortage  of  such  data,  an  empir¬ 
ical  approach  (based  on  radiative  transfer  simulations)  was 
applied  to  correct  contributions  from  Raman  scattering, 
similar  to  the  work  of  Loisel  and  Stramski  [2000]. 

[13]  Using  Hydrolight  [Mobley,  1995],  simulations  with 
and  without  Raman  scattering  were  carried  out  for  chloro¬ 
phyll  a  concentrations  equal  to  0.01  and  0.03  mg/W  (a(440) 
between  0.008  and  0.013  m  {),  and  for  solar  zenith  angles 
equal  to  10°  and  70°,  respectively.  From  these  simulations, 
the  ratio  of  R „  without  Raman  scattering  (R^no  Raman)  to 
/?„  with  Raman  scattering  (R^wah  Raman)  was  calculated  (rr  = 
^rs.no  Raman ^rs.wiih  Raman)-  ThcSC  rr  WCrC  3TOUnd  0.96,  0.95, 
0.89,  0.85,  and  0  84  for  410,  440,  490,  550,  and  670  nm, 
respectively,  and  were  quite  stable  (within  5%)  for  these 
simulations.  On  the  basis  of  these  analyses,  a  correction  of 
Raman  scattering  to  ScaWiFS/MODlS  measured  R was 
taken  as  rr  times  satellite  Rn,  and  this  product  was  then  used 


with  the  analytical  algorithm  for  further  data  processing. 
Without  such  a  correction,  the  derived  bbp  (via  analytical 
algorithms)  could  be  elevated  by  as  much  as  60%  for  such 
clear  waters. 

23.2.  Inversion  of  lOPs 

[u]  Previous  studies  [Gregg  et  al.,  2005;  McClain  et  al., 
2004;  Palter  et  al ,  2005,  Polovina  et  al.,  2008]  regarding 
the  oceanic  gyres  were  based  on  the  empirically  derived 
chlorophyll  a  concentration  from  the  observation  of  ocean 
color  [O Reilly  et  al.,  2000].  However,  such  empirical  algo¬ 
rithms  do  not  separate  the  optical  contributions  of  gelbstoff 
from  that  of  phytoplankton  pigments.  Consequently,  gelbstoff 
signals  may  be  '‘seen”  as  phytoplankton  [Carder  et  al.,  1 989]. 
To  delineate  the  phytoplankton  property  with  higher  fidelity, 
we  used  a  semianalytical  algorithm  [Lee  et  al.,  2002,  2007]  to 
derive  the  absorption  coefficient  (at  443  nm)  of  phytoplank¬ 
ton  (tfph)  and  that  of  the  combination  of  detntus  and  gelbstoff 
(adg),  respectively,  as  well  as  the  particle  backscattcnng 
coefficient  (/>b p).  Note  that  this  />bp  is  a  collective  index  rep¬ 
resenting  all  other  scattering  contributors  [Stramski  et  al., 
2001, 2004]  aside  from  water  molecules  (or  “pure". seawater) 

[is]  Because  molceular  scattering  makes  the  primary 
contnbution  to  remote  sensing  reflectance  in  the  blue-green 
wavelengths,  the  effects  of  the  molecular  scattering  phase 
function  need  to  be  well  accounted  for  [Lee  et  al.,  1994,  2004; 
Morel  and  Gentili,  1991;  Sathyendranath  and  Platt,  1997]. 
We  therefore  used  a  modified  analytical  function  for  R ^  [Lee 
et  al.,  2009]  based  on  the  model  of  Lee  et  al.  [2004]  that 
explicitly  separates  the  phase  function  effects  of  molecular 
and  particle  scattering 

+  (<»«>  + <" 

with  hi  =  a  +  />b,  and  bh  bhw  +  />bp.  Here  ft  represents  gen¬ 
erally  the  Sun-sensor  angular  geometry  for  the  remote  sens¬ 
ing  reflectance  and  taken  as  that  of  normalized  water-leaving 
radiance  in  this  study,  a  and  bh  arc  the  total  absorption  and 
backscattcring  coefficients,  respectively,  while  bhw  is  the 
backscattcring  coefficient  of  pure  seawater  (all  m  ').  Values 
of  the  model  parameters  (Go  (ft),  Gftft),  Go(ft),  and  G?(ft); 
sr  ])  for  various  Sun  angles  and  viewing  geometries  have 
been  developed  based  on  Hydrolight  simulations  [Lee  et  al., 
2009].  For  normalized  water  leaving  radiance  [Gordon, 
2005],  the  G  parameters  arc  0.0604,  0.0406,  0.0402,  and 
0.1310  sr  1  when  particle  scattering  takes  a  blended  phase 
function  [International  Ocean -Colour  Coordinating  Group 
(IOCCG),  2006]. 

[i6]  In  the  inversion  process,  the  backscattcring  and 
absorption  coefficients  of  pure  seawater  arc  based  on  the 
work  of  Zhang  et  al  [2009]  and  Pope  and  Fry '  [1997], 
respectively.  />bW(A)  is  calculated  with  a  salinity  value  of 
38.4%o  [Zhang  et  al.,  2009].  When  the  total  absorption 
coefficient  was  decomposed  into  the  contributions  of  phy¬ 
toplankton  and  dctritus/gclbstoff  [Lee  et  al.,  2002],  the  ratio 
of  tfpb(4!2)/flph(443)  =  0  742  was  used  for  such  low  chloro¬ 
phyll  waters  [Bricaud  et  al ,  2010],  and  the  ratio  of  ndg(412)/ 
fldg(443)  =  1.592  was  used,  which  was  calculated  with  a 
spectral  slope  of  0.01 5  nm  [Morel  et  al .,  2007;  Twardowski 
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Figure  3.  For  measurements  made  in  the  South  Pacific  gyre  (November  2004),  remote  sensing  inversions 
arc  compared  with  values  determined  independently  by  other  means,  (a)  Total  absorption,  (b)  phytoplankton 
pigment  absorption,  (c)  dctritus/gclbsto ff absorption,  and  (d)  particle  backscattcring  coefficient  at  550  nm 
(values  representing  a:  axis  were  estimate  based  upon  chlorophyll  a  concentration,  see  section  3.1).  Values 
for  jt  axis  were  based  on  Table  I  in  the  work  of  Morel  et  al.  [2007] 


et  al ,  2004,  Bricaud  et  ai ,  2010]  when  spectral  adg  is  mod¬ 
eled  as  an  exponential  decay  function  of  wavelength  [ Bricaud 
et  al l,  1981].  This  spectral  slope,  which  differs  from  that 
reported  by  Morel  et  al.  [2007]  because  of  differences  in 
wavelength  ranges,  is  based  on  measurements  in  NOMAD 
[Werdell  and  Bailey ,  2005]  (also  see  http://www.ioccg.org/ 
groups/So  ft  warc_OC  A/Q  A  A  v  5  pd  f). 

3.  Results  and  Discussion 

3.1.  In  Situ  Data 

[17]  We  used  data  published  by  Morel  et  ai  [2007]  to 
evaluate  the  inversion  performance.  Figure  3  compares 
derived  a( 420),  aph(420),  and  adg( 420)  with  those  presented 
in  Table  1  in  the  work  of  Morel  et  al.  [2007].  aph(420)  and 
adg(420)  that  represent  sample  measurements  were  converted 
from  the  measured  ap( 420)  and  ag(420)  [Morel  et  al. ,  2007], 
with  ad(420)  considered  10%  of  ap(420)  (aph(420)  then 
equals  ap(420)  -  ad(420)),  and  this  ad(420)  was  added  to 
tfg(420)  to  make  adg(420)).  The  10%  value  was  estimated 
based  on  previous  reports  [Siegel  et  al. ,  2002;  Morel  et  al , 
2007]  that  detritus  absorption  makes  a  small  portion  of  par¬ 
ticulate  absorption  for  such  extremely  clear  waters. 


[is]  As  described  in  detail  by  Morel  et  al.  [2007],  the 
optical  values  arc  extremely  low  for  this  oligotrophic  water  in 
November  (only  a  couple  of  months  after  the  winter  bloom 
[McClain  et  al ,  2004]),  with  a(420)  in  a  range  of  0.006- 
0.012  m  \  aph(420)  between  0.0017  and  0  0037  m  \  and 
tfdg(420)  between  0.0014  and  0.0036  m  1  These  extremely 
low  values,  all  from  sample  measurements,  arc  the  funda¬ 
mental  reasons  why  the  water  appears  “violet  blue”  [Morel 
et  al. ,  2007].  The  root- mean -square  differences  in  linear 
scale  between  /^-derived  and  published  absorption  coef¬ 
ficients  arc  -0.0010  m  \  0.0005  nT1,  and  0.0005  m  1  for 
<?(420),  aph(420),  and  adg(420),  respectively.  These  differ¬ 
ences  arc  equivalent  or  smaller  than  the  measurement  preci¬ 
sion  for  existing  instniments/mcthods  for  oceanic 
measurements  [ Twardowski  et  al. ,  2007].  On  a  relative  scale, 
the  averaged  percentage  differences  (average  of  I  measured 
/?re-dcrivcdl/measurcd)  for  the  three  optical  properties  arc 
-1  1%,  -18%,  and  -20%,  respectively.  If  ad(420)  is  consid¬ 
ered  as  5%  ofap(420),  the  average  percentage  differences  for 
aph(420)  and  tfdg(420)  become  -20%  and  -26%,  respectively. 
Considering  that  values  for  both  the  x  and  y  axes  are  asso¬ 
ciated  with  uncertainties  [Lee  et  al. ,  2010],  these  compari¬ 
sons  confirm  that  inherent  optical  properties  derived  from 
measured  R ^  arc  as  good  as  those  determined  by  direct 
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Figure  4.  Same  as  Figure  3,  but  for  comparison  of  Chi  a. 
The  OC4v4  algorithm  was  used  for  the  derivation  of  Chi  a 
from  reflectance. 

measurement  [Morel  et  al .,  2007],  at  least  to  the  measure¬ 
ments  made  in  November  2004.  Inversion  algorithms  can 
therefore  derive  highly  reliable  absorption  coefficients  from 
/?„  for  such  extremely  clear  waters,  as  has  been  recently 
shown  true  for  other  oligotrophic  and  mcsotrophic  environ¬ 
ments  [Gordon  et  al .,  2009].  This  robust  conclusion  is  par¬ 
ticularly  important  for  applying  such  algorithms  to  satellite 
data  to  obtain  long-term  observations  of  these  optical  prop¬ 
erties  of  such  waters,  although  additional  uncertainty  from 
satellite  measurements  will  be  introduced  [Antoine  et  al ., 
2008].  For  waters  in  the  open  ocean,  especially  gyres,  less 
uncertainty  (<I0%)  in  satellite -measured  R ^  is  expected 
because  of  the  stronger  w  ater-leaving  radiances  in  the  blue- 
green  wavelengths  and  the  easier  handling  of  atmosphere 
corrections  [Antoine  et  al .,  2008;  Gordon  and  Wang,  1994] 

[19]  The  ^-derived  particle  backscattenng  coefficient 
at  550  nm  (6bp(550))  is  approximately  0  0006  m" 1  for  the 
five  stations  (see  Figure  3d).  This  is  higher  (by  -0.0003  m  l) 
than  the  particle  backscattcring  coefficient  inferred  from  the 
measured  chlorophyll  a  concentrations  [Huot  et  al .,  2008; 
Twardowski  et  al ,  2007].  In  addition,  the  two  sets  of  bbp  are 
not  correlated.  Reasons  for  these  differences,  although  they 
arc  in  fact  within  the  measurement  precision  [Twardowski 
et  ai,  2007],  include  uncertainties  in  both  in  situ  measure¬ 
ment  and  remote  sensing  retrievals  as  well  as  other  factors 
that  arc  not  yet  included  in  the  inversion  process.  Such  factors 
include  effects  of  particle  phase  function  shape  [Mobley’  et  al ., 
2002],  scattering  from  bubbles  [Zhang  et  al. ,  1998],  and 
wave-focusing  cfTcets  [Zaneveld  et  al .,  2001].  Methods  to 
correct  these  secondary  effects  arc  still  under  development. 

[20]  For  the  same  data  set,  when  Chi  a  derived  from  R ^ 
with  the  OC4v4  algorithm  [O  Reilly  et  al. ,  2000]  is  compared 
with  Chi  a  from  water  samples  (Figure  4),  the  average  per¬ 
centage  difference  is  -37%,  which  is  a  greater  relative  dis¬ 
crepancy  than  that  of  aph(420)  and  adg(420)  (much  higher 
than  that  of  a(420)).  It  appears  that  most  of  the  difference 
came  from  a  systematic  bias  (underestimation  by  -0.0 1  mg/m3) 
of  remotely  sensed  Chi  a  (coefficient  of  determination  r2  = 
0.94).  This  may  not  be  surprise  as  the  coefficients  employed 
by  the  QC4v4  (or  OC3m)  algorithm  were  derived  empirically 
so  that  a  “best  fit"  was  achieved  for  a  “global”  data  set  that 


covers  various  regions  at  various  temporal  cycles  [O  Reilly 
et  al. ,  1998]. 

3.2.  Satellite  Data 

3.2.1.  Time  Series  of  Bio-optical  Properties 

[21]  With  the  above  validated  semianalytical  algorithm, 
monthly  <ar(443),  aph(443),  and  tfdg(443)  of  the  study  site 
were  derived  from  the  monthly  R ^  time  scries  measured  by 
ScaWiFS  and  MOD1S,  respectively.  The  particulate  back- 
scattering  coefficient,  bhp( 551),  was  also  derived  from  both 
measurements,  and  they  are  generally  in  a  range  of  0.0004 
0.0007  m  !.  However,  because  of  the  relatively  large  uncer¬ 
tainty  between  retrieved  and  measured  />^p  (Figure  3d),  we 
omitted  the  evaluation  of  bbp  temporal  variations  in  this  study, 
although  generally  an  inverse  correlation  was  found  between 
/^-derived  apb  and  /^-derived  bhp.  Figure  5  presents  the 
inverted  absorption  properties  between  January  1998  and 
December  2008,  with  thick  symbols  indicating  values 
between  July  2002  and  December  2007,  a  period  with  over¬ 
lapping  measurements  from  both  ScaWiFS  and  MOD1S.  For 
the  three  absorption  properties  ( a ,  aph,  and  adg;  wavelength 
dependency  is  suppressed  for  clarity)  derived  from  both 
ScaWiFS  and  MOD1S,  the  coefficients  of  determination  (r) 
between  ScaW'iFS  and  MOD1S  inversions  arc  0.93,  0.75,  and 
0.84,  respectively,  and  the  average  percentage  differences  arc 
3.6%,  12.5%,  and  16.4%.  These  differences  arc  significantly 
smaller  than  the  seasonal  variations  of  these  properties 
(a:  -45%,  apb :  -80%,  and  0dg:  -300%)  The  slightly  larger 
percentage  difference  for  <;dg  results  partially  from  the 
extremely  low  values  in  the  denominator  and  a  larger  dis¬ 
crepancy  of  /?re(412)  in  the  Austral  summer  of  2007  (sec 
Figure  2) 

[22]  For  the  Chi  a  values  derived  empirically  from 

SeaWiFS  and  MOD1S  (Figure  5d),  the  average  percentage 
difference  is  15.2%  (r2  0.81)  for  the  July  2002  December 

2007  period.  The  larger  percentage  difference  and  smaller 
r2  value  (compared  to  the  results  of  a(443))  suggest  a  slight 
discrepancy  introduced  by  the  two  empirical  algorithms 
adopted  for  each  sensor.  Nevertheless,  as  the  differences  from 
sensors  and  algorithms  are  much  smaller  than  the  seasonal 
variations,  these  results  indicate  that  the  bio  optical  proper¬ 
ties  obtained  from  the  two  sensors  arc  not  only  consistent  in 
their  temporal  patterns,  but  also  in  agreement  in  their  absolute 
magnitudes. 

[23]  Because  of  this  consistency,  the  following  analyses 
and  discussions  regarding  the  temporal  changes  of  these  bio- 
optieal  properties  arc  primarily  based  on  a  merged  time  senes 
from  ScaW'iFS  and  MOD1S  retrievals.  On  the  other  hand, 
because  SeaWiFS  measurements  covered  1998  to  2007, 
discussions  based  solely  on  ScaWiFS  measurements  arc  also 
provided.  The  data  merging  approach  is  similar  to  that  of 
Maritorena  and  Siegel  [2005]  [also  sec  JOCCG,  2007],  and 
the  1 1  year  merged  monthly  time  senes  is: 

'  P  =  P( ScaWiFS)  in  and  before  June  2002; 

P  =  (P(SeaWiFS)  -f  P(MOD!S))/2,  bclwccn  July  2002 

(2) 

and  December  2007; 

,  P  =  P(MODIS),  in  and  after  January  2008. 
where  P  represents  <arph,  <?dg,  and  Chi  a. 
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Figure  5.  Comparison  of  inverted  bio-  optical  properties  from  SeaWiFS  (solid  squares)  and  MODIS  (open 
circles)  R^. 


3.2.2.  Seasonal  Variations 

[24]  To  demonstrate  the  seasonal  dynamics  of  these 
bio-optieal  properties,  merged  monthly  properties  and  the 
monthly  time  series  of  SST  arc  presented  in  Figures  6a-6c. 
Clearly,  all  optical  properties  and  SST  show  well 
pronouneed  seasonality,  but  eaeh  has  its  own  timing  with 
respeet  to  seasonal  highs  and  lows,  as  highlighted  in 
Figures  6c  and  6d. 

[25]  The  phytoplankton  absorption,  tfph(443),  peaks 
(—0.0045  m  )  in  the  months  of  June  July  (Austral  winter) 
and  troughs  (—0.0025  m  ')  in  the  months  of  Dcccmbcr- 
January  (Austral  summer),  consistent  with  the  temporal  pat¬ 
terns  of  the  empirically  derived  Chi  a  (Figure  5d)  and  those 
from  earlier  studies  [McClain  et  ai.  2004;  Winn  et  ai ,  1995]. 
However,  the  aph  peaks  do  not  coincide  with  the  SST  troughs 
(August-September),  but  rather  occur  before  SST  reaches  the 
winter  minima  (Figure  6a).  This  relationship  agrees  with 
general  phytoplankton  dynamics  where  increases  in  biomass 
are  associated  with  enhanced  vertical  mixing,  which  delivers 
the  deepwater  nutrients  [McClain  et  al.,  2004],  as  June 
through  August  is  a  transition  period  with  SST  in  a  strong 
decline.  This  may  explain  a  low  coefTieient  of  determination 
(r2  =  0.09)  between  aph  and  SST  during  the  1998-2008  period 
(a  comparable  result  is  found  between  SeaWiFS  aph  and  SST 
for  the  1998-2007  period).  Also  note  that  the  aph  troughs 
follow  the  SST  minima  by  approximately  2  months  [McClain 
et  al .,  2004],  suggesting  a  combined  effect  of  nutrient 


depletion  and  reduced  solar  radiation  on  phytoplankton 
photosynthesis.  More  detailed  descriptions  on  seasonal  var¬ 
iations  of  phytoplankton  in  subtropical  gyres  ean  be  found 
in  the  work  of  McClain  et  al.  [2004] 

[26]  Different  from  the  relatively  mild  seasonal  variation  of 
tfph*  flag  shows  much  stronger  seasonality  (-0.0008  m  1  in 
summer  to  -0.003  m  1  in  winter.  Figure  6b)  at  this  location. 
Such  strong  seasonal  variations  of<;dg  will  in  turn  impact  the 
phytoplankton  dynamics  as  adg  directly  modulates  UV  radi¬ 
ation  in  the  upper  water  column  [Cullen  and  Neale%  1994; 
Tedetti  et  ai ,  2007].  Separately,  ddg  not  only  shows  sharp 
increases  and  decreases  around  the  Austral  winter,  the  sea¬ 
sonal  highs  and  lows  of  «dg  do  not  coincide  with  that  of  aph 
(Figure  6c),  consistent  with  the  observation  of  a  phase  delay 
found  in  North  Atlantic  waters  [Hu  et  al. ,  2006].  The  inde¬ 
pendent  adg  seasonality  can  be  further  manifested  with  the 

ratio  [Hu  et  ai ,  2006;  Jolliff  et  ai ,  2008],  where  a 
wide  range  of  seasonal  variation  ( — 0.3  -0.8)  is  found  for  this 
location  (see  Figure  6d).  As  discussed  previously,  the  tfdg/tfph 
ratio  is  not  constant  in  time  [Hu  et  al ,  2006;  Jolliff  et  ai , 
2008]  even  for  open  ocean  waters,  although  to  the  first- 
order  gelbstoff  absorption  was  found  to  covary  with  chloro¬ 
phyll  concentration  globally  [Morel,  2009]. 

[27]  Separately,  a  strong  relationship  was  found  between 
adg  and  SST  (r2  0.64,  not  shown;  it  is  0,63  for  SeaWiFS  adg) 
As  SST  is  directly  related  to  solar  radiation,  the  adg  seasonal 
variations  clearly  show  that  adg  at  this  location  is  regulated  by 
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Figure  6.  For  the  1 998  2008  period,  monthly  time  scries  of  (a)  aph  and  (b)  adg,  and  their  contrast  with  sea 
surface  temperature  (SST;  red  square  symbol  with  dotted  line)  (c)  The  phase  difference  between  aph  and 
adg.  (d)  Seasonal  variations  of  tfdg/aph  and  its  relationship  with  SST. 


phytoplankton  degradation  as  well  as  solar  photoblcach- 
ing  [Del  Vecchio  and  B lough,  2002;  Hu  et  al. ,  2006;  Nelson 
et  al.,  1 998;  Siegel  et  al,  2002;  Twardowski  and  Donaghay, 
2002]. 

3.2.3.  Interannual  Trends  From  Statistical  Analysis 
and  Impacts  of  Time  Frame 

[2s]  Following  the  approaches  in  earlier  studies  [McClain 
et  al,  2004;  Polovina  et  al,  2008],  an  intcrannual  trend  of 
relevant  properties  can  be  drawn  with  linear  statistical  anal¬ 
ysis.  To  minimize  the  effects  of  seasonal  modulation  on  the 
evaluation  of  intcrannual  variations,  the  intcrannual  trend 
was  analyzed  with  annual  mean  of  each  property.  Statistical 
results  from  such  analyses  arc  presented  in  Table  1,  with 
indicating  an  increasing  trend  and  indicating  a  decreasing 
trend.  It  is  found  that  there  is  a  decreasing  trend  for  aph  (r  = 
0.20)  for  the  1 998-2008  period  (r2  =  0  40  for  the  1 998-2007 
ScaWiFS  aph,  also  a  decreasing  trend),  but  an  increasing 
trend  for  adg  (r2  =  0.30)  for  the  same  period.  The  declining 
trend  of  aph  is  consistent  with  the  evaluations  of  Gregg  et  al 
[2005]  for  the  South  Pacific  gyre  based  on  empirically 
derived  chlorophyll  a  concentrations  (also  see  Table  1  and 
section  3.2.4).  This  apparent  declining  trend  is  to  a  large 
extent  because  of  the  higher  ap ^  (or  Chi  a)  values  in  1998, 
which  resulted  from  the  strong  1997-1998  El  Nino  event 
[Behrenfeld  et  al,  2001 ;  Kahru  and  Mitchell ,  2000,  Susanto 
and Marra,  2005].  Excluding  the  observations  of  1 998,  there 
is  no  apparent  trend  of  aph  for  the  1999-2008  period  (r2 
0.04).  The  increased  values  of  tfdg  arc  mainly  a  result  of  the 


higher  adg  value  in  2008  (resulting  from  an  abnormal  drop 
in  MODIS-measured  /?re( 412),  see  Figure  7).  Omitting  the 
adg  value  of  2008,  also  no  apparent  trend  was  found  for 
tfdg  (r2  "  0.06). 

[29]  Clearly,  sueh  “trends”  are  not  unequivocal.  Either  an 
increase  or  a  decrease  over  time  (trends,  and  with  a  different 
coefficient  of  determination)  eould  be  deduced  when  differ¬ 
ent  time  frames  were  used  for  the  statistical  analysis  (sec 
Table  I).  Note  that  sueh  different  trends  arc  not  sensor 
dependent  as  statistical  results  with  products  from  ScaWiFS 
alone  arc  nearly  the  same  (sec  Tabic  1,  except  that  aph  from 
ScaWiFS  shows  no  apparent  trend  for  the  1999  2006 
period.).  These  results  further  demonstrate  equivocal  inter¬ 
pretations  ofbiogcoehemical  responses  to  changes  of  climate 
when  using  sueh  short  records.  This  assertion  is  further  stated 
given  the  results  of  the  following  analyses,  which  are  focused 
on  the  merged  time  series. 

[30]  We  separated  the  monthly  time  scries  of  a  bio  optical 
property  into  annual  background  and  annual  intensity.  The 
annual  background  is  defined  as  the  average  low  value  during 
Austral  summer,  while  the  annual  intensity  (I)  is  defined  as 
the  difference  between  the  average  high  value  during  Austral 
winter  and  the  average  trough  during  previous  summer. 

I  =  Pcak(in  winier)  -  Trough(in  previous  summer).  (3) 

For  aph  and  Chi  a,  the  summer  and  winter  months  arc 
December-  February  and  June-August,  respectively.  Because 
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Table  1.  Trend,  Coefficient  of  Determination  (r2),  and  Significance 
Between  Time  and  aph,  fldg,  Chi  a ,  and  SST  for  Four  Separate  Time 
Frames'* 


Property 

1998  2003 
(P) 

1999  2006 
(P) 

2000  2005 
(P) 

1998  2008 

(P) 

apt>(443) 

0.09  (0.56) 
0.27  (0.29) 

+;  0.15  (0.34) 

;  0.00  (0.95) 

+;  0.67  (0.05) 
+;  0.35  (0,22) 

0.20  (0  17) 

ajg(443) 

0,24  (0.32) 

0  02  (0.80) 

+;  0.65  (0.02) 
+;  0.34  (0  13) 

+;  0.49  (0.12) 
+;  0.63  (0.06) 

+;  0,30  (0.08) 

Chi  ab 

0.22  (0.35) 
0.13  (0.49) 

+;  0.17  (0.31) 
+;  0.04  (  0.64) 

+;  0.67  (0,05) 
+;  0.55  (0.09) 

;  0  05  (0.52) 

SST  + 

0.04  (0.69) 

:  0.17  (0,31) 

;  0.63  (0.06) 

+,  0.00  (0.93) 

"In  the  table  +,  increase;  ,  decrease;  SST,  sea  surface  lcmpcraiurc.  The 
second  row  for  each  bio-optieal  property  indicalcs  evaluation  results  wilh 
products  solely  from  ScaWiFS. 

bChl  a  here  used  values  derived  from  the  empirical  algorithms 

tfdg  has  a  different  temporal  phase  from  that  of  aph,  the 
summer  and  winter  months  are  Fcbruary-April  and  July- 
September,  respectively. 

[31]  For  the  1998-2008  period,  the  annual  background 
absorption  of  ap h  is  ~0.0026  m  1  with  an  intensity  of 
''D.0013  m  1  (see  Figure  8a).  The  background  and  intensity 
of  tfph  also  showed  different  interannual  variations,  where  a 
more  pronounced  trend  emerged  for  the  background  of  ap h 
(a  decrease,  with  r2  =  0.24),  but  not  for  the  intensity  of  aph 
(an  increase,  with  r2  =  0. 1 1 ).  On  the  other  hand,  no  significant 
interannual  trend  was  found  for  the  background  of  adg  (sec 
Figure  8b),  although  the  adg  intensity  suggested  an  increasing 
trend.  The  interannual  variations  of  the  ap h  background  and 
tfph  intensity  suggest  that  even  if  there  is  a  decrease  of  phy¬ 
toplankton  at  this  location  in  the  1998-2008  period,  the 
decline  has  happened  during  the  “quiet”  period,  not  during 
the  period  when  temperature  decreases  and  vertical  mixing 
occurs.  Therefore,  the  dynamics  that  arc  responsible  for  the 
seasonal  increase  in  phytoplankton  pigments  has  apparently 
been  maintained  over  the  past  decade. 


[32]  These  observations  add  perplexity  to  the  interpretation 
regarding  declining  phytoplankton  in  oceanic  gyres.  Earlier 
studies  [Gregg  et  al ,  2005;  Polovina  et  al .,  2008]  showed  an 
intcrannual  decline  in  averaged  chlorophyll  a  concentrations 
for  surface  waters  in  the  South  Pacific  gyre,  and  speculated 
that  this  reduction  was  a  response  to  increased  water  tem¬ 
perature  (consequently  enhancing  stratification  and  reducing 
nutrient  supply  from  vertical  mixing).  However,  McClain 
et  al  [2004]  pointed  out  that  lor  this  subtropical  gyre  (and 
other  gyres  [ Palter  et  al. ,  2005;  Winn  et  a l ,  1 995]),  the  depth 
of  the  nutriclinc  is  significantly  deeper  than  the  mix  layer 
depth  (^250  versus  MOO  m)  so  the  supply  of  nutrients  might 
be  mainly  associated  with  surface  lateral  advection  (also  see 
discussions  in  the  work  of  Ayers  and  Lozier  [20 1 0]).  Also,  the 
temperature  time  series  did  not  show  a  trend  (r2  =  0.00)  at  this 
site  for  the  1998-2008  period  (see  Table  1).  In  such  a  sce¬ 
nario,  it  is  less  likely  that  the  reduction  of  phytoplankton 
pigments  in  the  surface  center  of  the  South  Pacific  gyre  results 
from  enhanced  stratification,  especially  when  the  reduction 
seems  most  apparent  in  the  quiet  period  of  each  year. 

[33]  Again,  it  is  necessary  to  emphasize  that  intcrannual 
trends  derived  from  statistical  analyses  are  sensitive  to  the 
time  windows  used.  For  instance,  if  we  focus  on  the  years  of 
1999-2006  (8  years)  or  200(F  2005  (6  years),  aph  actually 
revealed  an  increasing  trend  (see  Table  1 ),  further  suggesting 
limitations  to  infer  long  term  expectations  of  oceanic  phy¬ 
toplankton  from  short  observational  records  [Henson  et  al. , 
2009;  Martinez  et  al.y  2009]  Such  equivocal  conclusions 
might  be  caused  by  apparent  intcrannual  oscillations  (with  a 
cycle  of  ~6  years,  see  Figure  8)  in  the  background  and 
intensity  of  aph  (also  Chi  a%  section  3  2.4),  and  consequently 
both  decreases  and  increases  could  be  observed  when  short, 
but  different,  time  windows  were  used.  Dedicated  studies 
with  long-term  (multiple  decades),  consistent,  and  well 
calibrated  observations  arc  required  to  fully  account  for  the 
physical  forces  responsible  for  the  intcrannual  variations. 
Such  research  will  only  then  provide  the  fundamentals  to 
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Figure  7.  Monthly  time  senes  of  ratios  of  /?rs(443)//?Pi(551)  and  flrs(412)//?rs(443),  where  R ^  were  mea¬ 
sured  by  MOD1S.  For  this  study  site,  an  abnormal  drop  in  2008  Austral  winter  spring  (circled)  is  found 
for  the  /?rs(412)//?rs(443)  ratio,  which  indicates  larger  adg  for  similar  aph.  Further  study  is  required  to  confirm 
and  validate  this  drop. 
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Figure  8.  Interannual  variations  of  background  and  intensity  of  (a)  aph  and  (b)  adg. 


understand,  and  then  predict,  longer-term  secular  trends 
[Henson  et  aiy  2009]. 

3.2.4.  Comparison  With  Empirically  Derived  Chi  a 

[34]  Although  time  series  of  phytoplankton  absorption  can 
be  used  to  analyze  the  dynamics  of  phytoplankton,  tradi¬ 
tionally  such  evaluations  were  based  on  the  time  series  of 
empirically  derived  Chi  a  from  satellite-measured  ocean 
color.  Figure  9  presents  the  monthly  Chi  a  (derived  with 
standard  empirical  algorithms)  and  tfph  values  (both  are 
merged  results  of  ScaWiFS  and  MOD1S  products)  for  the 
study  site  over  the  1998-2008  period  Note  that  ap h  could  be 
converted  to  Chi  a  if  the  chlorophyll-specific  absorption 
coefficient  (a  *h)  is  known.  We  omitted  the  comparison  of  Chi 
a  here  because  for  the  same  Chi  a  value  aph  vanes  widely, 
especially  for  low  Chi  a  samples  (by  a  factor  of  ^3  for  low  to 
high  values  [Bricaud  et  al .,  1995]),  and  it  is  not  clear  yet 
how  to  best  predict  aph  value  remotely.  Instead,  the  time 
scries  of  Chi  a  and  aph  arc  compared  and  contrasted  here  In 
addition,  although  the  empirically  denved  Chi  a  is  on  aver¬ 
age  ~37%  lower  than  sample  measured  Chi  a  (section  3.1, 
Figure  4),  this  difference  is  assumed  to  have  no  impact  on 
analyses  of  temporal  variations  or  trends  as  it  represents 
mostly  a  bias  (slope  M.O  and  r2  ~  0.94  in  linear  regression) 
for  this  location 


[35]  Generally,  the  two  lime  series  show  nearly  iden¬ 
tical  seasonal  modulations  (except  October  December  2008 
where  a  higher  resulted  from  low  MOD1S  7^(4 12)  caused 
smaller  aph)  and  intcrannual  variations  (also  see  Table  1). 
This  observation  suggests  that  products  from  both  data- 
processing  approaches  arc  consistent  in  representing  phyto¬ 
plankton  signal  (and  then  intcrannual  trends)  at  this  location. 
This  consistency  results  from  the  facts  that  ( 1 )  Chi  a  derived 
with  the  OC4v4  (or  OC3m)  algorithm  (/?Pi(443)//?ni(55x)  is 
the  algorithm  input  for  oceanic  waters)  actually  represents 
the  total  absorption  at  443  nm  for  such  waters  (r  is  0.97 
between  QAA  derived  d(443)  and  empirically  derived  Chi  a , 
see  Figure  10)  and  (2)  the  absorption  of  CDOM  and  detritus 
makes  relatively  small  contributions  to  the  total  absorption 
at  443  nm  (Figure  bd). 

[36]  However,  there  are  differences  between  the  two  time 
senes.  First,  the  empirically  denved  Chi  a  shows  a  stronger 
seasonal  Chi  a  intensity  (see  Figure  1 1,  a  near  doubling  of 
Chi  a  during  the  winter  bloom  [McClain  et  al. ,  2004]).  The 
apparent  stronger  seasonal  bloom  is  partially  due  to  the  sys¬ 
tematic  bias  of  the  empirically  denved  Chi  a  (underestimation 
by  0.01  mg/m3,  see  Figure  4)  and,  consequently,  the  summer 
Chi  a  from  R „  could  be  too  low.  When  a  value  of  0.0 1  mg/m3 
is  added  to  the  time  series,  the  relative  seasonal  intensity  is 


'e 

a 


Figure  9.  Merged  monthly  time  series  ofChl  a  (blue  solid  squares)  and  aph  (open  circles),  with  aph  slightly 
peaking  ahead  of  Chi  a. 
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Figure  10.  Relationship  between  Quasi- Analytieal- 
Algorilhm -derived  tf(443)  and  empirically  derived  Chi  a  of 
the  study  site,  with  both  properties  retrieved  from  monthly 
Rrs  measured  by  SeaWiFS  and  MODIS. 

reduced.  In  addition,  because  of  the  impact  of  ODOM  and 
detntus,  the  strong  seasonal  elevation  of  adg  (see  Figure  5b)  is 
perceived  as  “growth"  of  phytoplankton  in  the  empirical 
algorithm.  A ftcr  remov  ing  the  effect  o f  a  dg,  the  aph  time  scries 
does  not  show  such  a  strong  winter  bloom  (Figures  5a  and  8). 
The  milder  seasonal  variation  can  also  be  deduced  when  is 
converted  to  Chi  a  with  an  a*ph(443)  value  of  0.1 1  mg/m2, 
where  Ohl  a  is  found  to  increase  by  just  ~70%  during  winter 
bloom.  This  a*Ph(443)  value  is  derived  with  ap  and  Chi  a 
values  in  Table  1  of  Morel  et  al.  [2007],  with  aph(443) 
approximated  as  ap(420).  This  a*ph(443)  value  is  also  con¬ 
sistent  with  the  average  of  tf*ph(443)  for  Chi  a  =  0.03  mg/m’ 
[Bricaud  et  al .,  1995].  Further,  the  intcrannual  variations  of 
background  and  intensity  derived  from  the  Chi  a  time  scries 
show  weaker  decline  trends  (sec  Figure  1 1). 

[37]  The  contrast  between  the  time  scries  of  aph  (derived 
scmianalytically)  and  Chi  a  (derived  empirically)  further 


demonstrates  that  if  an  accurate  description  (and  thus  better 
understanding  and  forecasting)  of  the  dynamics  of  phyto¬ 
plankton  is  required,  then  it  is  necessary  to  separate  phyto¬ 
plankton  from  gclbstoff  and  detntus  in  ocean  eolor  remote 
sensing  [Carder  et  al.%  1991]  even  for  such  extremely 
oligotrophie  waters. 

4.  Summary 

[38]  For  a  site  in  the  surface  center  of  the  South  Pacific  gyre 
(the  clearest  waters  in  the  global  ocean)  and  based  on  mea¬ 
surements  of  remote  sensing  reflectance  from  both  SeaWiFS 
and  MODIS  (Aqua),  we  produced  an  1 1  year  monthly  time 
scries  of  absorption  coefficients  of  phytoplankton  and  detri- 
tus/gelbstoff  and  a  time  scries  of  chlorophyll  a  concentra¬ 
tion  (from  standard  empirical  algorithms).  The  ocean  color 
inversion  process  for  the  optical  properties  is  composed  of 
steps  to  correct  the  contributions  from  Raman  scattering  and 
validated  with  in  situ  measurements  made  around  the  South 
Pacific  gyre.  The  time  senes  of  the  optical  properties  show 
clear  seasonalities  and,  more  importantly,  their  seasonalities 
follow  their  own  distinct  temporal  phases.  Such  results 
further  demonstrate  that  different  physical  and  biological 
processes  are  responsible  for  the  seasonal  and  intcrannual 
variations  in  the  optical  properties  of  this  hypcroligotrophic 
ocean  region. 

[39]  The  1 1  year  monthly  time  senes  suggest  a  decline  in 
phytoplankton  pigments  and  an  increase  of  dctntus/gclbstoff 
over  the  1998-2008  period,  although  statistically  the 
significances  arc  mild.  This  phytoplankton  pigment  trend  is 
consistent  with  that  depicted  by  empirically  derived  Chi  a% 
and  with  reports  by  McClain  et  al  [2004]  and  Gregg  et  al. 
[2005]  suggesting  that  changes  arc  happening  not  only  at 
the  fringe,  but  also  at  the  surface  center  of  the  South  Pacific 
gyre.  Previous  studies  [Gregg  et  ai ,  2005;  Polovina  et  ai. 
2008]  speculated  that  the  decrease  of  chlorophyll  in  the 
gyres  is  a  response  to  reduced  nutrient  supply  due  to  increased 
stratification.  However,  the  results  here  (especially  the  time 
scries  of  aph  annual  background  and  tfph  annual  intensity) 
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Figure  11.  As  in  Figure  8,  intcrannual  variations  in  the  background  and  intensity  of  Chi  a  (derived  from 
the  empirical  algorithms,  sec  text). 


11  of  13 


(09012 


LEE  ET  AL.:  BIO-OPTICAL  PROPERTIES  IN  A  GYRE 


C09012 


demonstrate  a  more  eomplex  picture.  This  is  because  (1)  the 
statistical  results  (or  interannual  trends)  are  sensitive  to  the 
time  windows  used  for  regression  analysis,  (2)  nutrient  sup¬ 
ply  is  less  likely  from  deepwater  [McClain  el  al. ,  2004],  and 
(3)  the  observed  reduction  of  phytoplankton  pigments  mainly 
happened  during  the  quiet  period  where  significant  stratifi¬ 
cation  is  already  evident.  Therefore,  in  addition  to  obtaining 
interannual  trends  of  biogeoehemieal  properties  [Anloine 
el  a! .,  2005;  Gregg  el  al .,  2005;  Polovina  el  al ,  2008],  it  is 
equally  or  more  important  to  identify  the  physical  forces 
and  quantify  their  relative  contributions  for  the  seasonal  and 
interannual  variations  [Ayers  and  Lozier ,  2010;  Marlinez 
el  al .,  2009;  Paher  el  al .,  2005;  Siegel  el  al .,  2008]. 

[40]  Because  of  the  relatively  isolated  location,  changes  in 
biogeoehemieal  properties  at  the  surfaec  center  of  the  gyres 
eould  provide  important  information  about  the  status  of  the 
oeeans  and  climate.  Consequently,  long-term,  consistent,  and 
reliable  measurements  of  such  properties  at  this  location  are 
an  important  part  of  the  monitoring  of  global  oeeans.  For  this 
objective,  semianalytieal  algorithms  that  separate  the  con¬ 
tributions  of  gelbstofT  and  phytoplankton  to  the  observed 
ocean  color  are  elearly  more  appropriate  than  simple  empir¬ 
ical  algorithms  (which  mainly  represent  ehanges  of  the  total 
absorption  coefficient),  especially  for  delineating  the  status 
and  dynamics  of  the  resident  phytoplankton  community. 
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